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Abstract 

We present a new method to address multipolar resonances and to control the scattered field of a spherical 
scatterer. This method is based on the engineering of the multipolar content of the incident beam. We 
propose experimentally feasible techniques to generate light beams which contain only a few multipolar 
modes. The techniques uses incident beams with a well defined component of the angular momentum 
and appropriate focusing with aplanatic lenses. The control of the multipolar content of light beams allow 
for the excitation of single Mie resonances and unprecedented control of the scattered field from spherical 
particles. 
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Resonances are ubiquitous in the Physical Sciences and Engineering. In a general sense, a reso- 
nance expresses that a system is efficiently coupled to an external excitation within the interval of a 
certain parameter. When properly controlled, this efficient coupling can be used to measure differ- 
ent phenomena with high sensitivity. For instance, the control of nuclear magnetic resonances in 
nuclei has enabled a great number of applications in biology and medicine yj]. Unfortunately, the 
lack of control of these resonant phenomena can derive in catastrophic results, one of the most fa- 
mous being the collapse of the Takoma bridge in the 1940's M\. Recently, resonances in plasmonic 
systems have been gaining a lot of interest [3] and some techniques to control them have been pro- 
posed [|4j-l6|]. Most resonances can be approximately modelled with the Lorentzian distribution. 
However, in 1961 a new fundamental kind of resonance was discovered by Fano [7]. Although it 
was firstly regarded as a quantum mechanical effect, the Fano resonance is very present in optics 
andplasmonics. In particular, one of the most important analytical solutions of Maxwell equations 
in inhomogeneous media, i.e. Mie Theory, gives rise to Fano resonances, too [8D. The Mie Theory 
isl or the more modern Generalized Lorenz-Mie Theory (GLMT) Jlp| describe the interaction be- 
tween an incident electromagnetic field propagating in a lossless, homogeneous, isotropic medium 
and a homogeneous isotropic sphere. The number of applications of this theory is immense. The 



reader is referred to Gouesbet et al. for an extensive review 111 ill 



The most appropriate electromagnetic modes for describing the physical phenomena associated 
with the GLMT are the multipolar modes: {A^, A^}, where (m) and (e) stand for magnetic 
and electric multipole. These modes are eigenvectors of the total angular momentum (AM) oper- 
ator J 2 and one of its projections such as J z with respective values j and m z Ill2l ll3n. This set of 
modes exploit the spherical symmetries of the kind of problems studied by the GLMT. As a conse- 
quence, the electromagnetic eigenmodes of the scattered field from a spherical object are precisely 
the multipolar modes. Experimentally, it is very challenging to prepare and detect a single mul- 
tipole with an arbitrary value of both j and m z . This is due to the spherical and non-paraxial 
character of these modes. In this paper we propose a technique to control multipolar resonances of 
dielectric spheres by engineering the multipolar content of experimentally realizable light beams. 
In our approach, we use focused cylindrically symmetric beams. These light beams possess a well 
defined value of J z . In 1992 Allen and co-workers showed that in the paraxial approximation one 
could find a set modes which had a well defined value of the AM fl^j . A particular set of paraxial 
modes with this property are the Laguerre-Gaussian modes (LG q j, with q the radial index and 
I the azimuthal one). This important finding opened up a whole new field, where the ability to 



independently control the AM and the polarization components of paraxial beams has allowed in- 
numerable applications in quantum optics, microscopy, biology, optical trapping and astrophysics, 
just to mention a few Jl5ll . The total AM of a vectorial field is composed of an orbital part (L) and 
a spin part (S), i.e. J = L + S. Nevertheless, it is important to note that these two components 
cannot ordinarily be separated. When both of these operators are applied to a Maxwell field, the 
result from that operation is not a Maxwell field [16, §16] [17]. 

This paper presents two main findings. First, we show how to experimentally control the mul- 
tipolar content of cylindrically symmetric modes. We advance that by controlling a component 
of the AM and the focusing of the incoming beam we can engineer fields which contain just few 
predetermined multipolar modes. Secondly, we demonstrate that by properly engineering these 
modes, we can also control the scattered field as well as address single resonances of an arbitrary, 
homogeneous and isotropic sphere in a lossless, homogeneous and isotropic medium. The second 
idea falls under the GLMT and is highly relevant to applications where control of the scattered 
field from a spherical object is desired. This scattering problem is solved by decomposing all the 
fields involved in the problem (incident E 1 , scattered E sca , and absorbed E abs ) into multipoles and 
then use the boundary conditions. The three fields of the problem have the same formal expres- 
sion. However, the multipolar decomposition of the fields E sca and E abs have all their complex 
amplitudes modulated by a Mie coefficient. These Mie coefficients are obtained by applying the 
Maxwell boundary conditions and depend on the AM as well as the parity. Thus, the only techni- 
cality of the problem is the projection of the incident field into the basis of multipoles. 

The cylindrical symmetry of the incident beam simplifies the problem since the other two fields 
are required to possess a well defined value of J z , too. In this way, we can solve the problem using a 
mixture of analytical and numerical techniques. More importantly, by exploiting the symmetries of 
the system we can get a further insight on the problem only using considerations of conservation of 
AM. Previous works considered similar problems where spherical metallic particles were excited 



with Laguerre Gaussian beams [18]. There, the use of purely numerical techniques allowed the 
authors to distinguish certain properties of the scattered field. Our approach is different and more 
general. We extend the previous results to excitation with any cylindrically symmetric field and 
show a clear way of controlling the scattering from spherical particles. The results presented in 
this paper also allow for a better understanding of recent experiments of scattering of silica spheres 
with Laguerre-Gaussian beams Jl^l . as our formulation also includes the paraxial case. 

The experimental set-up which we propose is as follows. First, we prepare a monochromatic 



paraxial beam with a well defined J z controlling its polarization and spatial properties separately. 
This can be typically achieved with either spatial light modulators or holograms and a set of 
waveplates and polarizers. After this step, we will focus the resulting beam with a lens which 
we model as aplanatic I20I . Chap. 3]. The resulting focused beam interacts with the sphere and 
we analyze the scattering field. From now on, we will express our results in terms of the vector 
potential in the Coulomb gauge. All electric fields can be recovered using the typical formulas, i.e. 



E = ikA 1120. Hence, the expression of the paraxial beam previously mentioned is A'(p, <f>,z = 
Z\) = Fi(p) exp(il(f)) exp(ikzi)e p , where {p, 0, z} are the spatial cylindrical coordinates, Z\ is the 
transverse plane right before the lens, k is the modulus of the wavevector, Fi is an arbitrary function 
representing the complex radial amplitude and e p = (x + i py) / \/2 is a unit polarization vector 
representing right (p = +1) or left (p = —1) circular polarization. Here, we are disregarding terms 
of the order 1 / (kd) 2 , where d is the typical size of our light beam, which in our case is of the order 
of the focusing lens diameter. Any paraxial beam can be decomposed as a superposition of this set 



of modes |21|]. Note that a paraxial beam with circular polarization has a well defined helicity A, 
with a very good approximation [22]. As mentioned before, then the paraxial beam impinges on 
the back of the microscope objective and it is focused. Microscope objective manufacturers make 
the objectives so that they behave approximately as aplanatic lenses. In addition to this, they coat 



the lenses in a way such that the transmission coefficients for the s and p polarized waves [23, 



Chap. 13.1] are equal, i.e. t s {6) = t p {6). Therefore, specifying the aplanatic model described in 
I21I Chap. 3] for this case, and taking into account that this model preserves the helicity and the 
AM of the incident beam 11711 ■ we can write the field at the focus of the lens as 
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A foc (r) = C / d(f) k / sin 9 k d9 k \/n cos 9 k x 
Jo Jo 

F(f sin 9 k )exp(i(l + p)(j) k )M(9 k , (f) k )e p exp(ikr) (1) 

where C = f exp(—ikf)/27i; f is the focal of the lens; 9m = arcsin(NA), with NA the numerical 
aperture of the lens; k = k(sin 9 k cos 4> k , sin 9 k sin <j) k , cos 9 k ); n is the index of refraction of the 
lens; r = (pcos(</>), psin(0), z) is the position vector; and M(# fc , (f) k ) is the rotation matrix applied 
to the polarization vector so that the transversality of the field is fulfilled. Note, that field (OQ) fulfills 
the exact Maxwell equations. 

We will now describe how these beams will allow us to control the scattered field off a sphere. 
As previously mentioned, this problem is solved once the incident field (A foc in our case) is de- 
composed into multipolar modes. This decomposition can be made using the method put forward 
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FIG. 1. Multipolar decomposition (|Cj/ p | 2 ) for different cases. The insets represent the intensity plots of the 
modes used for each simulation. The red color indicates NA=0.25, and the blue NA=0.9. The multipolar 
decomposition of a) LGrj,o an d b) LGo,3 is presented. Note that Cji p can be described with very few 
multipoles when we use a high NA microscope objective to focus the beam. 
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j=\l+p\ 



A (m) 4- 7nA (e) 
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Cji p = C / sin 9 k d9 k d j {l v) {9 k )Fi (/ sin 6 k ) y/n cos 9 k 
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As previously stated, the lens does not break the cylindrical symmetry of the beam, meaning that 
after passing through the lens, the beam still carries the same AM component J z = m z . Thus, 
one can see in the decomposition © that m z is the same for all multipoles and corresponds to 
m z = l+p, which bounds the value of j > \l+p\. The coefficients of the decomposition Cji p only 
depend on a one dimensional integral which averages the radial complex amplitude of the input 
beam (F{) over the reduced rotation matrix d J mp (9), which can be found in Il25ll . Finally, note that 
the beam still has a well defined helicity with value p. Indeed, if we apply the helicity operator 
fl26|] A = (Vx)/fc to A foc in eq. © we obtain AA foc = pA foc . This is a consequence of the 
following two relations |12I1 : 
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Fig. 1 exemplifies one of the main results of this paper. It shows the multipolar decomposi- 
tion of different focused beams. Two different input beams were used: a Gaussian beam and a 
Laguerre-Gaussian with optical charge I = 3, both right circularly polarized. The width of both 
beams was chosen so that the entrance pupil was filled. We consider two aplanatic lenses with 
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FIG. 2. Suppression of background in the scattering efficiency (Q sca ). Input beams are a) plane wave and 
b) LGo,i8- In a) no aplanatic lens is used, whereas in b) an aplanatic lens with a NA=0.9 is used. The rest of 
parameters are kept constant for both plots, i.e. r = 1.3 fim, n r = ^Je r ^i r = 1.5, p = 1. The insets indicate 
the typical profile of the beam used to excite the sphere. The scattering efficiency is represented with a blue 
continuous line. The Mie coefficients 620 and 022 are plotted with a red and a green line respectively. In a) 
we have indicated with a dashed line the position of two particular resonances for these modes. Note that 
the ordinate axis in a) and b) are different. 

NA's equal to 0.25 and 0.9. Our results show that the higher the NA, the narrower the distribution 
of multipolar modes is. The other parameter to control the distribution of multipoles is the J z 
of the incident beam. As mentioned earlier, the distribution of multipolar amplitudes is zero for 
values of j < \l + p\. This allows us to control the multipolar content of an optical field. 

The decomposition © and the use of the results of Mie scattering are the needed elements to 
calculate the scattered and interior fields, 



A sca = _ £ ij(2j + l)y 2 C jlp 

j=\l+p\ 
oo 

A abs = Y, ^(2i + l) 1/2 C, 7p 
j=\l+ P \ 



h-X {m) +7W7-A (6) 
°3 A j(l+p) + l P a 3 A j(l+p) 
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where {aj, bj, Cj, dj] are the Mie coefficients defined in II27L Chap. 4]. These coefficients only 
depend on the size parameter of the problem (x = 2nr/\), and both the relative permeability 
(fx r ) and permittivity (e r ) of the sphere with respect to the surrounding medium, r is the radius 
of the particle. This completes the solution of the problem, as we can now easily compute the 
value of the fields in all regions. In particular, we can calculate the scattering efficiency, which 
is a dimensionless measure of the scattered power by the sphere in the far field 11271 Chap. 3]. 
Using the incident beam given by eq. ©, the following expression is obtained for the scattering 
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efficiency: 



Q™= £ *^|cWW + M a ) (5) 



i=|f+p| 



In comparison, the expression of the scattering efficiency for an incident plane wave 112 8l Chap. 
10], has all | Cj-^j, | 2 equal to 1 and the sum starts at j = 1. This means that in conventional Mie 
Theory, all the Mie coefficients (and consequently all the multipolar modes) have the same relative 
weight. In our case, the use of an engineered input beam allow us to control the relative weight of 
the different multipolar modes (given by aj and bj) with the use of the multipolar content of the 
incident field (given by Cji p ). Then, the a tightly focused incident beam will allow us to select a 
few relevant Mie scattering coefficients. On the other hand, the use of beams with a well defined 
value of J z allows us to exclude the first \m\ multipolar modes from the scattering cross section. 

In order to show these effects, in Fig. 2(a) we plot the conventional Mie scattering efficiency 
when the incident field is a plane wave, as well as the values of a pair of high order Mie coefficients 
(620 and 022)- In this case, the cross section is dominated by the low order multipoles. The higher 
order resonances appear as little ripples on top of the large oscillation. This is the so-called ripple 
structure. However, if an engineered cylindrically symmetric beam is used, the lowest \m\ modes 
can be removed and thus a particular Mie coefficient of interest can be excited. This is shown in 
Fig 2(b), where the scattering efficiency associated to a LG ,is beam is depicted, as well as the 
same two Mie coefficients plotted in Fig. 2(a). 

One interesting possibility which raises from the ability to freely excite multipolar resonances 
is that one can address very sharp resonances. The Mie scattering coefficients oscillate from 
to 1 in a very exotic fashion. This behaviour is given by spherical Bessel functions. Every 
resonance of a Mie coefficient is given by a zero of a transcendental equation involving Bessel 
functions of order j. These equations have an infinite number of zeros, giving raise to an infinite 
number of resonances for a given order j. In general, it is very difficult to compare the shape 
of different resonances. Nevertheless, the first peaks of the Mie coefficients do follow a general 
trend. The higher the order of the Mie coefficient, the sharper the resonance. This implies that 
the high multipolar resonances of the sphere have smaller linewidths and consequently very large 
Q factors. This also means that this sphere will be much more sensitive to the incident radiation. 
This fact is used in order to excite whispering gallery modes in large spheres, reaching very high Q 
factors for incredibly high values of j [29|]. In theory, Q could reach values as high as 10 400 . Yet, 



due to radiative losses of the material, the maximum values of Q experimentally achievable are of 
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FIG. 3. Effect of a) the input beam and b) the NA of the lens on the excitation of a single resonance in the 
scattering efficiency (Q sca ). In a), three different beams are used: a plane wave and two LG's beams with 
I = 4, 8 focused with a lens with NA=0.25. In b) we use a LGo,i8 (same as Fig. 2(b)) but we focus it with 
different NA=0.1, 0.3, 0.5 and 0.9. The beams always fill the entrance pupil of the lens. Note that in a) the 
ordinate axis is linear and in b) logarithmic. Also, in a) the efficiency produced by each of the LG beams is 
multiplied by a different factor so that all the curves are comparable in the same plot. 



the order of Q = 10 9 It30tl . Although this feature should also be present in the scattering of light 
from spheres, it is typically hidden due to the background produced by the lower order modes. 

As shown in Fig. 3(a), our technique allows to excite a single multipolar resonance with free- 
space radiation. We fix our attention to the little ripple created by the Mie coefficient 6 2 o m the 
plane wave scattering (also present in Fig. 2(a)). The figure shows how to reduce the background 
created by the excitation of the lower modes. This has a power drawback though. Increasing the 
AM of the incident mode reduces the absolute value of the scattering cross section if the NA of 
the lens is preserved. However, if we increase the NA, the retrieved scattered power increases by 
many orders of magnitude and the shape of the resonance is maintained. This is shown in Fig. 
3(b). Thus, high NA microscope objectives in conjunction with high AM modes can be used to 
experimentally address single multipolar resonances. 

In Fig. 4 we show the absorbed and the total electric field (scattered plus incident) surrounding 
the sphere for two different incident beams. Two sets of four plots corresponding to two different 
incident beams (LG ,o and LG 0i is) at the same combination of {r = 1.3 fj,m,n r = 1.5} are pre- 
sented. Both incident fields have a well defined helicity A + (or p = 1) in the focus of the lens. 
Thus, the resonant behaviour of the scattered field will be specially clear when the opposite com- 
ponent of the helicity A is examined, since the contribution due to the incident field is 0. Indeed, 
if we compare the resonant excitation shown in Fig. 4(f) with the off -resonant behaviour of Fig. 
4(h), it can be seen that the intensity differ in two orders of magnitude. Note that as it was shown 
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FIG. 4. Projection into the helicity basis of the intensity of the absorbed and total field from a sphere with 
parameters {r = 1.3 /j,m,n r = 1.5}. The helicity of the incident beam is always p = 1. The intensity is 
plotted in logarithmic scale, where corresponds to the maximum of intensity in a). Images a), b), e) and 
f) have been simulated for a A reS onance = 503 nm, whereas a 2 nm shift has been introduced in the off- 
resonance wavelength for the others, i.e. A ff. reS onance = 505 nm. The simulations for the LGo,o beam have 
been executed with a NA=0.25 lens. On the contrary, a lens with a NA=0.9 has been used for the excitation 
with a LGo,is- In both cases, the entrance pupil of the lens was filled. 



in Fig. 2, the resonance depicted here is driven by the 620 Mie coefficient, hence the single excited 
multipolar mode is A 2 ™^ 19 . Another consequence of the single excitation of a multipolar mode in 
Fig. 4(e-h) is the fact that in Fig. 4(e-f) the two helicity components have the same power. This is 
due to the fact that the we are exciting a single magnetic multipole A^ 19 and that the multipoles 
considered in © are a superposition of multipoles with a well defined helicity II 1711 . Finally, it can 
be observed that unlike Fig. 4(e-h), where the field can be described with a single mode, many 
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multipoles are needed to describe the interferences seen in Fig. 4(a-d). In fact, the pattern-like 
observed in Fig 4(a-h) has drawn a lot of attention since 2004, where it was first characterized as 



photonic nanojetJ31|]. Since then, nanojets have been extensively studied and applied in many 
different fields [13211. 

We have demonstrated how to address multipolar resonances and how to control the scattered 
field of a spherical object by engineering the incident light. This technique allows to excite single 
multipolar resonances and also, to choose resonances with very high Q factors even for spheres 
with sizes of the order of the wavelength of the incident radiation. Besides the fundamental control 
of material particles that this technique grants, we foresee applications in the fields of cytometry 
and dark field microscopy, where the control of the AM of the illuminating beams has not been 
explored until now. 

This work was funded by the Australian Research Council Discovery Project DPI 10103697. 
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